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1. Introduction 

In the ground state of electron systems, it has been shown that exchange (X) due to the 
Pauli 'repulsion' and correlation (C) due to the Coulomb repulsion suppress particle- 
number fluctuations and consequently reduce the energy This energy reduction 
provides most of the 'glue' that binds atoms together to form molecules and solids @]. 
Particle-number fluctuations mean that the particle number in a domain (which may be a 
muffin-tin sphere, a Wigner-Seitz cell, a Bader basin [f|, a Daudel loge ||, a bond region 
between atoms in a molecule, etc.) fluctuates due to zero temperature quantum motion 
with a certain probability. Fulde Q takes C2H2 as an example for such fluctuations. 
The number of valence electrons in a sphere containing a C atom fluctuates around 
its average value ~ 3.9. Comparison of Hartree-Fock (HF) calculations for C2H2 with 
calculations which include correlation shows that configurations with large deviations 
from the average valence electron number (e.g., with and 1 or 7 and 8 electrons) 
are strongly suppressed due to correlation. A similar fluctuation-correlation analysis is 
performed in Ref. for several dimers and in Ref. || for the uniform electron gas in one, 
two, and three dimensions (ID, 2D, 3D). These calculations for the above mentioned 
narrowing of the particle-number distribution need the pair density (PD) nifi,^) and 
this narrowing is used to derive from the PD a quantum-kinematical measure for the 
correlation strength |T[. Correlation and its strength is furthermore characterized by 
the small-separation (or on-top) behavior of the PD. The spherically averaged on-top 
curvature of the spin-parallel PD may serve as a local correlation measure |7|] and from 
the topological analysis of the intracule PD a short-range correlation strength is defined 
||. In addition to these PD based quantities the concept of a correlation 'entropy' has 
been developed for Fermi systems |9Hl2f (in Ref. the term Jaynes entropy is used). 



It is based on the correlation induced non-idempotency of the correlated one-particle 
density matrix (1PDM) / y(r;r f ). All these correlation measures intend to make the 



qualitative terms 'weak and strong correlation' quantitatively precise [ 13| . Note that 
strong correlation means extreme narrowing of the particle-number distribution which 
is usually described as electron localization. 

In the present paper, we apply the above mentioned fluctuation-correlation analysis 



to the exactly solvable Calogero-Sutherland (CS) model Uli| . The CS model is a model 
of long-range-interacting spinless particles in ID and has been solved exactly by means 
of the (asymptotic) Bethe-Ansatz technique f!4|, |15||. The solution is valid for both 



fermionic and bosonic particle symmetry. Here we will mostly concentrate on the Fermi 
systems. Furthermore, the model can be shown to be the universal quantum model 
underlying the dynamical interpretation of random matrix theory |f[5|, IT||. This latter 
connection has been used to also compute several correlation functions exactly at three 
special values of the interaction strength, among them the 1PDM and the PD flTJJ . Thus 



although the information is restricted to the ID case, the model nevertheless is ideally 
suited for testing the fluctuation-correlation measures discussed above. 

Correlated 1PDM and correlated PD need correlated many-body wave functions 
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(beyond the HF caricature), which in quantum chemistry [T8|, [19j are traditionally 
obtained from configuration interaction (CI), coupled cluster (CC), M0ller-Plesset, 
quantum Monte Carlo calculations or recently from the contracted-Schrodinger-equation 
method |20| or the incremental method ||21|| . All these procedures involve certain 
approximations or have restricted applicability. So the existence of non-trivial exactly 
solvable models which can provide 1PDM and PD is of much interest for the mentioned 
fluctuation and correlation analysis. 



We shall only consider the ground state properties of the CS model |T4|. The 



interaction is pairwise inversely proportional to the square of the distance = \xi — x 3 , 



of two particles with interaction strength /z, i.e., p/xfy The interaction strength 



/i > —1/4 is occasionally parametrized as fi = A(A — 1) with a parameter A 



1/2 + a/1/4 + > 1/2. We shall mostly use the parameterization v = a/1/4 + /i, 
such that fi = v 2 — 1/4 and A = 1/2 + v. In the thermodynamic limit we assume 
constant density p[x) = n, so the CS-ground state has only two parameters, v and n. 
The interaction has the peculiarity not to possess a natural length. Therefore it is 
a model showing critical behavior, which can be discussed in terms of universality classes 
and their conformal anomalies p2]-p5|. This beauty of the 1/xf,- interaction shows up 
also in the analytical Bethe-Ansatz solutions fli"4|, p6H30| and the explicit knowledge of 
the correlated many-body wave functions |T4] , |31| . From the Bethe-Ansatz technique the 
complete energy spectrum and in particular the ground state energy per particle as a 



function of the interaction strength parameter v is available [|14[] We show that its kinetic 
and interaction 'components' can be deduced with the help of the Hellmann-Feynman 
theorem [52]. Surprisingly, when the interaction strength parameter v approaches its 
limiting value 0, both the kinetic and the interaction energy diverge in such a way 
that they compensate each other leaving the total energy finite. As we outline in the 
following, these divergences result from the peculiar behavior of the 1PDM and the PD 
for v — > and are related to the "fall-into-the-origin" already mentioned in Ref. [|33| . 

For v = 0, 1/2, and 3/2 — corresponding to /i = —1/4,0 and 2 or A = 1/2, 1, and 2 
- it has been shown ]14|] that the square of the ground state wave function is intimately 
related to the eigenvalue distribution of random matrices of the orthogonal ensemble, 
the unitary ensemble, and the symplectic ensemble, respectively. Using this connection, 
Sutherland had shown how to construct the 1PDM and the PD using integral relations 
of random matrix theory. The resulting formulas reduce the problem, say for the 1PDM, 
from the evaluation of a high-dimensional integral to the computation of a determinant 
of a matrix [15]. From the 1PDM j(x—x'), the momentum distributions n K for the three 
special values of v follow via Fourier transform. Due to correlation the latter quantities 
are non-idempotent. They determine the mentioned correlation entropy per particle 
s = — n K lnn K / n K . Knowledge of the PD n(x 12 ) allows |T3], [T3[ [|§] us to 



calculate the fluctuation ANx of the particle-number around its mean value Nx = nX 
in any piece (domain) X of the a;-axis. Comparing this variance of the particle-number 
distribution Px{N) for the cases 'no correlation' [y = 1/2 or HF approximation) and 
'correlation' shows the above mentioned narrowing for repulsion [y > 1/2) in a smaller 
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ANx- Contrarily, for attraction [v < 1/2) a broadening with a larger ANx appears. 

In Section 0, we introduce the CS model, define the kinematical quantities used 
throughout the text, and present the Hellmann-Feynman theorem. Section |3| is devoted 
to the thermodynamic limit. In Section [|, after presenting the HF approximation, we 
discuss first qualitatively and then analytically the influences of the CS interaction on 
1PDM and PD. In particular, we show that the above mentioned divergences in kinetic 
and potential energies are caused by a peculiar behavior of the PD n(xij) for small inter- 
particle separations Xij <C k^ 1 and of the momentum distribution n K for large momenta 
> fcp or k e k/kp 3> 1. In Section || we then apply the mentioned fluctuation- 
correlation measures to the CS model. Section || is devoted to details of the numerics 
and in Section |7] we discuss extensions of our approach to impenetrable Bosons and 
lattice gases. We conclude in Section || with a discussion of our results. 



2. The system and its ground state 

2.1. Hamiltonian, energies, and quantum kinematical quantities 
The Hamiltonian of the CS model is H = T + V with 

AT N n 2 _ 1 

f = H- 2 ti > * = +£^> (2.i) 

i i i<j V 

with pf = —d 2 /dx 2 , and N equal to the number of particles. We assume the 
system to be confined to the length L by an external potential v ext (x), e.g., a box 
or harmonic oscillator potential. In the following we alternatively assume periodic 
boundary conditions with v ext (x) = and a density in the k space described by 



LAk/{2n) = 1 [ I4"[ , P5 |. The average particle density is n = N/L. Furthermore, it follows 
from dimensional reasons that all energies for the Hamiltonian ( |2.1| ) are proportional to 
n 2 in agreement with the virial theorem and all lengths are measured in units of 1/n 
(thus kp ~ n) 0. 

We denote the ground state energy and its kinetic and potential 'components' by 
En = (H), T/v = (T), and Vn = (V), respectively. Then the corresponding energies 
per particle are e^ = E N /N, tjy = T N /N, = V^/N with e^ = t^ + v^. From the 
antisymmetric (or symmetric) ground state wave function follow by contractions [14], lq 



the 1PDM 7at(x;x') and the PD ^(21,22) normalized as J dx^N{x;x) = N and 
J dx\ J dx2 n^(xi,X2) = N(N — 1), respectively. The PD describes the XC hole, 
vanishing for zero separation and approaching the Hartree product Pn(xi)pn(x2) for 
large separations. The cumulant PD u7 A r(xi,x 2 ) = Pn(xi)pn(x 2 ) — n N (xi,x 2 ) is via 

dx\ \ dx 2 wn(xi,X2) = N (2.2) 



size-extensively normalized. Furthermore with the abbreviation y = k^x, kp = irn fT4" 
and with the dimensionless functions /W^i; 2/2) hermitian, gjy(yi, ^2) non-negative, and 
h N {yi,y 2 ) = /jv(yi;yi)/jv(z/2;y2) - ^(2/1,2/2), we can write for the 1PDM j N (xi;x 2 ) = 
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n /jv(?/i; 2/2)) f° r the PD n N (x\,x 2 ) = n 2 g^ivi, y 2 ), and for the cumulant PD we have 
Wn(xi,x 2 ) = n 2 h N (yi,y 2 ). The dimensionless cumulant PD is thus h N = 1 — g N and 
normalized as 



1 

N 



dyi dy 2 



7T 7T 



(2.3) 



which follows from Equation (|2.2j ). With these dimensionless 1PDM and PD and with 
the Fermi energy e F = fcf/2 the energies t^ and vn are given by 



1 

N 



dyi 

7T 



dy{ 



/at (y 1:2/2) 



V2=yi 



and 



9N{yi,y2)— cf- 

7T 2/12 



f2.4a) 



(2.46) 



Therefore t N /e F , fjv/ep and e^/e-p are functions of /x and iV. The latter dependence 
disappears for the thermodynamic limit as shown in Section |3|. 



2.2. The Hellmann-Feynman theorem 

If e^ is known as a function of /i, then t^ and Ujv can be obtained from the Hellmann- 
Feynman theorem [3l| without knowing the quantum-kinematical quantities /at(?/i; 2/2) 
and Pat (2/1 5 2/2)- This theorem says 



9/i \ d/jj / 
which for (|2.1|) gives 



;2.5) 



v N = /j,^-, t N = ( 1 - fi-^- ) e N (2.6) 



V d/j 



d d . ... . 

—t N = -Utt ~ v n ■ (2.7) 



and also 

djj, r 9/i \/i 

Thus — with Equation ( |2.1| ) in mind — the Hellmann-Feynman relation ( |2.5| ) for the 
l/x 2 j model establishes an integral relation between the dimensionless 1PDM /jy on the 
l.h.s. and the dimensionless PD on the r.h.s. of Equation 



3. Thermodynamic limit 

We wish to study the thermodynamic limit with iV — > 00 and L — > 00 such that 
n = N/L = const. The resulting extended system has only two parameters, the 
interaction strength parameter v and the Fermi wave number k F . So i/ep, v/ep, and e/e F 
become functions of v only. The thermodynamic limit makes furthermore the 1PDM 
and the PD to depend only on k F Xi 2 = k F (xi — x 2 ) = y 1 — y 2 = y 12 . The dimensionless 
functions fx, gN, and then take the forms f(yi 2 ), g(yu), and h(yi 2 ) = 1 — g{yu), 
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respectively, with /(O) = 1 and g(0) = or equivalently h(0) = 1. These functions have 
v as the only parameter. 

The eigenfunctions (or natural orbitals) of the 1PDM 7(2:12) = 71/(2/12) become 
simply plane waves y>\{x) = e lkx /L, such that 

POO 

7(^12) = n / dK n K cos K\y 12 \ = n f(y 12 ), (3.1) 

where n K is the momentum distribution and k = k/kp. 

For v = 1/2 (ideal spinless ID Fermi gas) the Pauli principle leads in the reciprocal 
space to the Fermi ice block n° K = 9(1 — \k\) and in the direct space to the ideal X 
hole g°(y) = 1 — [f°(y)] with the dimensionless 1PDM f°(y) = (sin y)/y following from 
Equation (|3.1|) . The energy per particle is e = ep/3 = /cf/6. 

In general, with 7(0) = n, the momentum distribution n K is normalized as 
n K = N or 



00 



dK n K — 1 . (3.2) 
The kinetic energy per particle is according to Equation ( |2 . 4 oj) 



K 2 



t = 6 dK n K —e . (3.3) 
Jo * 

n K is a function of \k\ and v, so t/e^ is a function of v only with t = e$ for v = 1/2. 
The corresponding expressions for the PD g(y) are according to Equation (|2T 



2 r*M.h{y) = l, h{y) = l-g{y) (3.4) 
Jo 71 

and for the interaction energy per particle according to Equation ( 2.4lf) 
r°° dy V + | 

g(y) is a function of y and z/, so f/eo is a function of z/ only. With t and t> follows the 
integral relation 

00 

, k 2 dn K f dy v 2 + \ dq(y) 



as a consequence of the Hellmann-Feynman theorem expressed in Equation (|2.7|) . 
Correlation via v 7^ 1/2 deforms the X hole and the Fermi ice block as shown in Figs. 
H] and § in such a way that Equation (|3.6| ) is maintained. 



4. Hartree-Fock approximation and correlation beyond it 

4-1. Hartree-Fock (HF) approximation 

The simplest approximation for the quantities n K , g(y), t, and v is obtained from the HF 
approach. The momentum distribution in Equation (|3.3| ) and the PD in Equation (|3.5|) 
are to be replaced by their 'ideal' expressions n° K and g°(y), respectively. Consequently, 
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we find £hf = and whf = 2/ie and thus chf = (1 + 2//) e , as shown in Figure [I]. Here 
the identity ( |A.5| ) has been used. The total HF energy eHF also obeys the Hellmann- 
Feynman theorem ([2.5|) and the virial theorem. 



4-2. Qualitative discussion of correlation 

Due to correlation the true ground state energy per particle, e, is below the HF 
energy enF and the true ground state wave function $(• • •) is no longer a single Slater 
determinant. Note that the definition of the term 'correlation' needs a reasonable 
reference state, which is $hf(" • •) in our case. So, correlation causes a negative 
correlation energy e corr = e — chf < 0, namely through redistributions of g°(y) and 
n° which are shown in Figs. [| and ^ and described in the following. 

As we show in Figure 0, correlation modifies the X hole of the unperturbed PD. 
Especially the correlation induced changes for small y are of interest, because the 
interaction 1/y 2 is there largest. The on-top behavior of the uncorrelated X hole 
(y = 1/2 or HF) is described by g°(y) = y 2 /3 + . . .. In its correlated counterpart 
with a //-dependent exponent and //-dependent coefficients (see |Appcndix ri| ) 

g(y) = Ay a (l + a x y + a 2 y 2 + •••), 

a =l + 2v, (4.1) 

correlation for v ^ 1/2 shows up in a ^ 2 and A ^ 1/3. More precisely, repulsive particle 
interaction [v > 1/2) supports the Pauli 'repulsion', so the X hole is broadened (through 
increasing a and decreasing A), but attractive particle interaction iv < 1/2) fights 
against (or competes with) the Pauli 'repulsion', so the X hole is narrowed (through 
decreasing a and increasing A) as shown in Figure and Table [l]. This X hole narrowing 
(for v < 1/2) or broadening (for v > 1/2) makes 

9 r*M = 1 + -L %2 for (4.2) 
lo n y 2 2v 2 

The equation follows from Equation ( |3.5| ) together with the Hellmann-Feynman theorem 

(p.5|). Thus v <i>hf = 2/ieo for v ^ 1/2 as shown in Figure [I]. 

As shown in Figure || correlation thaws the Fermi ice block #(1 — |k|). 

Mathematically, the momentum distribution n K becomes non-idempotent, physically 

this means: Correlation excites particles for k > 1 and holes for k < 1. This increases the 

kinetic energy independent whether the interaction is attractive [y < 1/2) or repulsive 

(y > 1/2): t > tjjF as can be seen in Figure 0. We note that n K has no discontinuity 

at \k\ = 1. Its value is 1/2 and near k — 1 it follows a power law as is typical for all 

Luttinger liquids with their z-p = |36|,[37 . 



We model the melting of the Fermi ice block analytically with the continuous 
function 

n K = \ + B(l - Kf [I + &r(l - ^ + b 2 (l - ^ + ...} 

for < k < 1, (4.3a) 
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n K 



B{K - if [1 + b+(K - if + bt{K - l) W + -..} 



for 1 < « < 2, 
'1 + ^ + ^ + 



for 2 < k < oo 



(4.36) 



(4.3c) 



with [15 



= 7 



1 (1 - 2i/) 



(4.4) 

4 l + 2i/ V ; 

and 7 = 3 + 2z/ ( |Appendix B| ). The exponents (3, 7 and the coefficients B, C, and 6*, q 
depend on v. The condition < (3 < 1 ensures an infinite slope of n K at k = 1. This n K 
has to obey the normalization (|3.2|) and the condition 



3 / c£k Tl K K 
Iq 



2u 



> 1 



(4.5) 



which follows from Equation Q3.3|) together with the Hellmann-Feynman theorem (|2.5|) . 
For v — 1/2 (or HF) it is /3 = 0, 5(1 + £\ bf) = |, and (7 = 0. The correlation induced 
melting for v ^ 1/2 shows up in /3 > 0, B{1 + £\ bf) < I and C > 0. 



4-3. Results of the exact solution 



With the help of the Bethe-Ansatz technique one obtains e = A cq |14],[15 



6 cLS cL 

function of the interaction strength parameter A shows no special behavior for A > 1/2, 
but as a function of the interaction strength \i = A(A — 1), 

I - • I e , v=\/- + tt (4.6) 



2 y «7 V 4 r 

the non- analytical behavior for \i — > — 1/4 is incorporated in the variable v. 

Equation (|4.6|) yields with the Hellmann-Feynman theorem (|2.6| ) the kinetic energy 
per particle, 



2v 



e . 



(4.7) 



Equation (|4.6|) yields with Equation (|2.6|) also the interaction energy per particle 



M 1 + 2^) CO- 

From Figure [l] we see that both t and v diverge for v— >0, while e remains finite. 
Equations (fO) and (O) lead to 



(4.t 



dy g{y) 



(4.9) 



/o t r 

as another integral relation between the momentum distribution n K and the 
dimensionless PD g(y) in addition to Equation (|3.6j ). These distribution functions have 
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to change with v in such a way that these relations ( |3.6| ) and ( |4.9| ) are obeyed together 
with the normalization conditions ( p.2| ) and ( |3.4j ). 

The PD n[x\2) = n 2 g(y) with y = kpX\2 is known analytically for the values v = 0, 
1/2, and 3/2 For v = it is (with the notation of [Appendix""^ 

g{y) = 1 -(^)\ Si{v) ^_iA.^i, (4.io) 

\ V J dy V %dy y 

for v = 1/2 (ideal Fermi gas) it is 

M-l- (4.11) 



y 



and for v = 3/2 it is 



9to) . 1 _(^) +Si(2 , ) _|_^. (4 , 2) 

The corresponding dimensionless cumulant PDs h(y) = 1 — g(y) are given in Table ^| 
together with their Fourier transforms 

POO 

h(q) = 2 dycos(qy)h(y). (4.13) 
Jo 

They have via S(q) = 1 — h(q)/n a simple relation to the static structure factor (or 
van Hove correlation function) S(q) = (p q pu/N, which describes the correlation of 
density-density fluctuations. p q = exp(— igxj) is the Fourier transform of the density 
operator p(x) = £\ S(x — Xi). The three PDs g(y) are shown in Figure || and the three 
structure factors S(q) in Figure f|. 

For v = 1/2 the weak oscillations of g(y) and the (first-order) kink of S(q) arise 
from the Fermi momentum distribution n K with its sharp discontinuity zp = 1 at k — 1. 
With increasing repulsion the oscillations of g(y) are enhanced, what is displayed in the 
reciprocal space by the peak of S(q) at q = 2 (and a 3rd-order kink at q — 4). This 
peak is customary in ID quantum liquids, see, e.g., the spin-correlation functions of 
Ref. f39 |. The first maximum of g(y) runs through a certain trajectory from (it, 1) to 
(2.99, 1.24). Whereas repulsion enhances the Friedel oscillations of g(y) and the kink 
structure of S(q), they diminish with increasing attraction: for v = both g(y) and 
S(q) approach the value 1 smoothly (non-oscillatory) from below. The discontinuity for 
v — of the second derivative of S(q) at q = 2 replaces the kink seen for other values 
of v. For the on-top behavior of g(y) in terms of g{0), g'{0), g"{0) the following holds: 
It is g(0) = 0, according to the Pauli principle, </(0) = n/6 for v = 0, but for v > 0, 
and it is g"(0) = for u = 0, infinite for < v < 1/2, but 2/3 for u = 1/2, and for 
v > 1/2. We remark that S(q) and g(y) can be computed also for other values of v in 
addition to the three special values used here [|34"| , |4"5|| . 

With the identities ( |A.2|) - ([A.4|) the normalization condition ( |3.4|) is fulfilled. From 
Equations ( |4.10| )- (|4.12| ) follow the on-top coefficients of Equation ( |4.1| ); they are shown 
in Table [I]. Note that the last two terms of Equation ( |4.10| ) do not contribute to the 
normalization because of Equation ( |A.3j ) and that the last term causes the odd on-top 
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coefficients of Table [I] and also the linear behavior for small y. Its oscillations are exactly 
canceled by the combined oscillations of the second and the third term. Simultaneously, 
these terms ensure the correct normalization. 

The PD ( [4. 11| ) for v — > 1/2 plugged into Equation ( [2.4 £| ) yields with the identity 
( |A.5| ) the same as results from Equation Q4.8| ), which follows from the total energy per 
particle, Equation ( [4.6| ), and the Hellmann-Feynman theorem ( |2.6| ), namely v = 2/^eo- 
Similarly the PD (|4.12| ) for is — > 3/2 plugged into Equation (|3.5|) yields with the identities 
( |A.5| ) and (|A.6|) the same as results from Equation (|4~8| )), namely v = 4/ieo/3. 

For v = a divergence appears, because from the PD ( |4.10| ) follows an on-top 
behavior, which is linear in y as shown in Figure |2| and Table p]. This linear behavior 
results from the last term of Equation (|4.10| ), which does not influence the normalization 
( |3.4|) , but it makes the interaction energy v ~ J °° dy g{y)/y 2 to diverge logarithmically 
in agreement with the divergence of v — > —eo/8is for y^O as displayed in Figure [l]. 

The divergence of the interaction energy is accompanied and compensated by the 
corresponding divergence of the kinetic energy t — > eo/8v. This indicates a special 
asymptotic behavior of the momentum distribution n K for namely Equation ([4.3c]) 
with 7^3. For 7 > 3 the integral J °° dn n K n 2 is convergent, but with 7^3 for is^O it 
diverges logarithmically, whereas the normalization integral ( |3.2j) remains convergent. 
The counterpart to this asymptotic behavior of n K for k — > 00 is the on-top behavior 
of the PD for y — > as shown in Figure ^| and Table [I] with a smooth transition of 
the coefficient A in Equation (|4.1|) from 7r/6 via 1/3, to 16/135 for is = 0, 1/2 and 
3/2, respectively. With quadratic interpolation of the coefficients shown in Table [I] as 
functions of is, one may continuously switch the on-top behavior of the PD g(y) between 
its forms at is — and 3/2. For the PD exponent a = 1 + 1v we refer to |Appendix P| , 
where also the momentum-distribution exponent is conjectured as 7 = 3 + 1v. 

These divergences of the kinetic and the interaction energies indicate that for 
attractive particle interaction /i/x?- with \l— > — 1/4 the system becomes unstable (no 
ground state with finite kinetic and potential energies). We remark that it was shown in 
Ref. [33| that the singular particle interaction — |/i|/|ri2| 2 makes already two particles to 
fall together ( "fall-into-the-origin" ) for > 1/4 (ground state with E — > —00) and for 
\fx\ < 1/4 there are only scattering states with E > (no bound states with E < 0) fl4]| . 

For v = the exact solution of the CS model yields the momentum-distribution 
data. In Section || we will give the details of the necessary numerical calculation. 
The coefficients of Equation ([4.3 dj ) are fitted to the n K values for n = ... 1 and the 
coefficients of Equation ( 4 . 3 6| ) are fitted to the n K values for k — 1 . . . 2. The coefficients 
are chosen to make also n K at k = 2 continuous and smooth. Finally b% is fine tuned 
to make the normalization equal to 1 according to ( |3.2| ). The results are shown in 
Figure [5] and the values of the coefficients are given in Table ^. The case v = 3/2 is 
similarly treated only with the difference that also the kinetic energy t = |eo can be 
used for the fine tuning in addition to the normalization condition ( |3.2| ). The results 
are shown also in Table [| Here b% and 64 have been used for fine tuning (which yields 
the normalization 0.997 and t = 1.34 e , instead of 1 and 4e /3). 
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5. Fluctuation-correlation analysis 

5.1. Quantities following from the pair density 

Particle-number fluctuations: Following Fulde [Q one may ask to what extent 
correlation influences particle-number fluctuations ANx in a domain X, i.e., a certain 
interval of the x axis, where in the average there are Nx = nX particles. These 
fluctuations are measured quantitatively by |],|3],[l3| 

\2 i r X r X 



Nx nX J J 



dx2 w(xu) 



U Y d Vl f dy^, (5.1! 
Y Jo Jo 71 



with Y = kpX = nnX. Following the Appendix A of Ref. || the 2D integral ( |5.1|) is 
reduced to a ID integral with the help of the Fourier transform (|4.13| ), namely 

(ANx) 2 , 2 [°°dq l-cos(qY) h(q) 
N X ~ Wo tt q- tt ■ {b - 2} 

The results are shown in Figure |5], where also the case v — > oo ('strict' or 'perfect' 
correlation) is displayed. Traces of the oscillations as v — > oo are already visible for 
v = 3/2. 

With h(y) = 1 — g(y) and with the expansion of g(y) according to Equation ( |4.1| ) 
- see also the text after Equation ([4.12 ) — the small-A expansion of Equation ( |5.1| ) is 



(ANx) 2 _ 1 , A „ v , J , VN2 , , 



, dinX + d 2 (nX) 2 + d z {nX) c 

Nx 

+ d 4 (nX) 4 + d 5 (nX) 5 + . . . . (5.3) 



The slope d\ at X = does not depend on the interaction strength parameter 
v as shown in Table f| because of g(0) = and h(0) = 1 not depending on u; 
but the coefficients of the next terms do such that the particle-number fluctuations 
are suppressed due to repulsive particle interaction, but enhanced due to attractive 
particle interaction: correlation makes the particle- number distribution Px(N) more 
narrow for repulsion (y > 1/2) and more broad for attraction (y < 1/2). We 
remark, that fluctuation enhancement (induced by attractive interaction) generally 
may support /cause clusterings (e.g., paramagnons prior the para-to-ferromagnetic phase 
transition). In our case this tendency shows up in the sudden "fall-into-the-origin" at 
v = 0. If one considers with X — 1/n a Wigner-Seitz 'sphere' (with 'radius' A/2 and 
N x = 1), then 

ElH = 1 _ *H xM =<^M! (5.4) 
U X-(l/2) " ' N x 1 ' 

is a reasonable correlation measure based on particle-number fluctuations as we show 
in Figure 0. 
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On-top behavior: The exponent a and the coefficients A, a, of Equation (|4.1|) describe 
the short-range or dynamical correlation, i.e., the small-separation behavior of g(y), 
see Table [1]. Cioslowski's correlation cage || is in our case simply the inter-particle- 
separation range y = . . . y max with y max being that separation where the PD g(y) has 
its first maximum g max = g(y max ). For v = 0,1/2,3/2 the corresponding values are 
?/ma X = oo, 7T, 2.99 and g max = 1, 1, 1-24 ||. One may ask to what extent the correlation 
cage contributes to the interaction energy and define 



< 1 



Wl/2)' ' W fo^dyg^/y* ^ (5 ' 5) 
as an energetic correlation measure with V^ age (0) = 1; the expression simplifies when 
using (|4.2|) . Both Ei and E 2 vanish for v = 1/2 as shown in Figure || 

5.2. Quantities following from the momentum distribution 

Critical exponent: The critical or correlation exponent j3 of Equation (|4.4j) can be 
computed from conformal field theory ||15|,|38|. It describes (together with the coefficient 
B) the behavior of n K near k = 1 according to Equations (|4.3o|) and ( [4.3 fl) . For the 
three special values v = 0, 1/2 and 3/2, this gives 1/4, 0, and 1/4, respectively, as shown 
in Figure ^. The exponent 7 describes the decay of the correlation tail. 



Non-idempotency and correlation entropy: The qth order non-idempotency is [|T^ 
c(q) = 1 — J °° dn {n K ) q . The derivative of c(q) at q = 1 is s = c'(l) or 



s{y) 



dn n K In n K > 



(5.6) 



to be referred to as correlation entropy [[T|,[T3|]. It has been plotted in Figure 



Correlation tail properties: The relative number of particles (or holes) in the 

(5.7) 



corresponding correlation tail is [12, 13, JC 

iVtail(f) = 



dn n K 



dn (1 — n K ) < 1. 



The contribution of the correlation tail to s is p~3|| 



S ta n(v) = - I dn n K \nn K < s(u). 



(5.J 



(5.9) 



In addition to these quantum-kinematic measures one may use |T§] 
. , f, 00 dtv n K n 2 

T te il(^) = J -foc~. — < 1 

as another energetic measure with T tai i(0) = 1. Also these correlation measures vanish 
for v — 1/2 as shown in Figure |6|. 
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5.3. The correlation energy 
For e corr = e — em? follows 



V 2 



e corr = - ( v - - 1 e . (5.10) 

Kinetic and interaction energy contribute t corr = — j^e COTT and v cor r = (l + ^7) e corr , 
respectively, to e cor r- 

5.^. Comparison of the correlation measures 

When comparing the computed correlation measures in Figs. ||| it turns out that for 
small \v — 1/2 1 the PD based measures Ex .2 of Eqs. (|5.4j) and ( |5.5|) are proportional to 
v — 1/2 (which is —e' COTT /(2eo)), whereas the n K based measures ( ^.6|) -( pl9|) as well as 



( |4.4|) behave like iy — 1/2) 2 (which is — e corr /e ). So the latter ones are not so sensitive 
as the first ones. With s(v) = 0.5828|e corr /e | + . . . the Collins' conjecture |e corr | ~ s 
is confirmed at least for weak interaction. In this limit also iVtaii, Stan and Ttaii are 
mutually proportional and their derivatives are proportional to Si and E2. 

We remark that the quantities x{ v ), K age (^), N tan (v), StaiiM, and ^tau(^) are 
reference free, i.e., they are defined without reference to the non- interacting case v — 1/2 
- which in our case is simultaneously equivalent to the Hartree-Fock approximation. 
Reference quantities appear in Ei j2 with x(l/2) and V CSuge (l/2) and in s(u) with 
s(l/2) = 0. Whereas this observation is important for quantum chemistry — as stressed 
by J. Cioslowski M — whenever multi configuration appears, it is less important in our 
case which is well described by single configuration. 

6. Numerical determination of 1PDM and n(«) for the CS model 

Results from the theory of random matrices enable the calculation of various correlation 
functions for the CS model |T^ , p^| ,p8H30|]. In particular, the 1PDM can be expressed in 



terms of a determinant of an appropriate matrix Fpq |i~4| , |i~5|| . The size of this matrix is 
specified by the number of particles N to be (N — l) 2 for v = 1/2 and 3/2 and (iV — 1) 2 /4 
for v — 0. Each element of contains simple trigonometric ID {y — 1/2 and 3/2) or 
2D [y = 0) integrals. 

For some cases, most notably v — 1/2, the resulting determinant can be computed 



analytically and corresponding expressions have been given in Ref. [15|. For the other 



cases, we have evaluated the determinant numerically fI5 |, using a subdivision of the 
system volume (periodicity length) according to L/L$ = 42,402, and 402 for v = 0, 
1/2, and 3/2, respectively. The particle number, odd due to periodicity of the wave 
function ||15|| , varied from AT = 1 to 401, corresponding to a variation in density n from 
nearly to nearly 1. Taking the Fourier transform, we next compute the momentum 
distribution n K for all densities. In the inset of Figure |^, we show results for one of the 
three special v values. 
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Next, we apply the definitions of correlation measures and correlation energies as 
given in Sections |3], |4], and |5] and study their density dependence. In Figure [7| we 
show results for the entropy s and in Figure |8| for the various energies as the density 
is varied. As all energies scale with n 2 , these measures should be density independent 
when normalized with respect to e$. However, we do in fact see a pronounced density 
dependence for n > 0.5/L and also for n < 0.05/L . This latter density dependence is 
simply due to the small particle numbers, thus a small size of Fpq and consequently a 
limited resolution when computing the 1PDM at fixed L/Lq. The density dependence 
at large n values is more intricate to explain. The computation of the 1PDM by the 
connection with random matrix theory works for the periodic model p5| . Thus there 



exists a Brillouin zone and the tail of n K for \k\ outside this Brillouin zone is folded 
back into it. The tail of n K thus tends to be dominated by this effect for large n values 
as shown in the inset of Figure [3[ However, knowing that the correlation measures 
must be independent of density in the thermodynamic limit, we deduce their values by 
restricting us to these density regions where the independence holds. Then we apply 



the fit according to Equation ( |4.2| ) as explained in Section |4.2| . In Figure |6| we indicate 
by error bars the small variation in the correlation entropy when using instead of n K as 
in Figure | the n K as in the inset. Similarly, the corresponding variations in iV tail , S ta n, 
and T ta ii are within the symbol sizes. 

7. Extension to impenetrable Bosons and lattice gases 

As mentioned in the introduction, the CS model is also solvable for bosonic particle 
symmetry. The bosonic wave functions have to obey an additional boundary condition, 
namely they have to vanish for inter-particle separations Xij — > such that the 
resulting system consists of impenetrable or hard-core particles |14| with additional [ijx 2 ^ 



interaction. Both PD and 1PDM may be calculated as before. The PD is independent of 



statistics |[14|| , thus the fermionic exchange hole agrees with the bosonic impenetrability 
hole and all quantities computed before based on the PD are the same in the bosonic 
and the fermionic case. For the 1PDM this is different, the momentum distribution of 
Bosons is quite different from the fermionic n K as shown in Figure |9]. However, energetic 
quantities and correlation measures based upon those are nevertheless independent of 
the statistics and should thus be the same for Bosons and Fermions. In Figure |8| we show 
that this is indeed the case. Thus besides the density independence we have another 
criteria that allows us to extract the correct values of the correlation measures from 
these plots. We note that the abovementioned unwanted density dependence is also 
present in the bosonic n K and visible in Figs. and Also present is the aliasing effect 
as shown in the inset of Figure 0. 



In Refs. W\ ITS], it had been shown to be useful to restrict the family of wave 



functions of the CS model for both bosonic and fermionic symmetry to a lattice such that 
the coordinates are integers Xj = 1,2, . . .L |4l|-443|. Only the normalization constants 



of the wave functions change and the 1PDM can be computed much as before fl4j] , 
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replacing the integrals in Fpq by appropriate sums Furthermore, the structure 

factor S(q) is known exactly and therefore also the PD [15|. The resulting lattice gas 
has a particle- hole symmetry and thus we need to consider n < 1/2L only. However, 
the density N/L now enters all expressions in a non-trivial way and the very useful 
density independence of the continuum model for the quantities considered here is no 
longer applicable. Nevertheless, the continuum model corresponds to the low-density 
limit of the discrete model. In Figure |7|, we show that this is indeed the case for, e.g., 
the correlation entropy. 

8. Discussion and Conclusions 



Both the PD based and the n K based correlation measures (44 ) and (|5.4|) - (|5.9[) vanish for 



v = 1/2 (no interaction). But the first ones are more sensitive because of S 12 ~ v — 1/2 
near to the no-interaction point as shown in Figure |] (c), while the second ones are 
~ {y — 1/2) 2 like e corr of Equation (|5.10|) as shown in Figs. | (a) and (b) and therefore 



cannot distinguish between attractive and repulsive interactions. In ID the PD based 
measures (5.4) and ( |5.5D are identical for fermionic and (hard core) bosonic particles. 



The n K based measures ([14]) and ( |5.6| )- (|5^ ) do not apply for bosonic particles, they 
are designed for fermionic particles only and appropriate bosonic variants still have to 
be defined. 

For repulsive particle interaction results well-known from other extended many- 
body systems, such as enhancement of the Friedel oscillations with maxima/minima 
trajectories (Figure g), humps/peaks of the static structure factor developing from its 
non- interacting kink (Figure |j), suppression of particle-number fluctuations (Figure |5]), 
are confirmed again. However, for attractive interactions, we have found in the 
present work that particle-number fluctuations are contrarily enhanced and that this 
is accompanied by a smoothening of the PD (the oscillations disappear) and of the 
static structure factor (the kink disappears) as well as by the appearance of a linear 
on-top behavior of the PD. The latter behavior results in a diverging interaction energy 
in the strong attraction limit although the total energy remains finite. In momentum 
space the Fermi ice block thaws for both cases and correlation tails develop. In the 
strong attraction limit the correlation tail becomes so long ranged that the kinetic 
energy diverges, thereby exactly compensating the divergence of the interaction energy. 
We have shown that these divergences can be derived from the exactly known energy as 
a function of the interaction strength with the help of the Hellmann-Feynman theorem 
( p.5| ). This theorem allows to calculate t{v) and v (u) separately from e[y) and gives 
- in addition to their normalizations (13. 21) and (13.41) — exact relations for n K and the 



PD as shown in Equations (|3.3| ) and ( |3.5| ) together with Equations ( |4.7|) and 
respectively. Equations ( |3.6| ) and ( |4.9| ) give additional exact integral relations between 
n K and the PD. We expect that similar relations can be derived for correlation function 



results [f4*4|, |45|j based on the theory of Jack symmetric polynomials for values v = p/q 



with p and q relatively prime positive integers |54|. Thus an extension of our work to 
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these v is possible, albeit necessating a different approach for the numerically stable 
evaluation of the products of integrals. 

In summary, we have analyzed particle-number fluctuations and studied measures 
for the correlation strength based on the pair density and on the momentum distribution 
of the ID quantum system of interacting particles. We made extensive use of 

the available exact solution and applied the Hellmann-Feynman theorem to the CS 
model. Our results show that further work is called for in order to make the qualitative 
terms 'weak and strong correlation' quantitatively precise presumably with a variety of 
quantities instead of a single index |13| . 
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Table 1. On-top exponent and coefficients of the PD according to Equation (4.1). 
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Table 2. Dimensionless cumulant PD h(y) and the structure factor S(q) used for 

') as in Equations (f 
S(q) = 1 - %)/tt 



the computation of ANx and £i(z/) as in Equations (5.2) and (5.4). 
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Table 3. Coefficients as in Equation (4.3 a) - Equation ( 4.3t| ) calculated from the 
numerically determined momentum distribution n K for v = and 3/2 (at n = 1/2Lq). 
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Table 4. Coefficients of the smaii-X expansion of ^ A ^^ as in Equation ( |5.3| ). 
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Appendix A. Certain integrals 

The following identities are valid with Si(x) = J* dy[sm(y) / y]: 
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Equations (|A.2| ) - (|A.4|) determine the normalization of the PD's ( 4.10Q - (|4. 12 



Equations (|A.5|) and ( |A.6|) determine the interaction energy v for the HF approximation 
and for u = 3/2. For the fluctuation analysis with Equations ( |5.1|) and (|5.2| ) 
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Appendix B. Kimball like theorems for n(x\2) and n K 

The small separation or on-top behavior of the PD n{x\2) is derived here similarly 
as Kato [[|7| and Kimball [48 found the cusp relation dg(kpr)/dr\ r=0 = g(0)/a-B [or 
g'(0) = ar s g(0),a = (4/97T) 1 / 3 ] for the pair correlation of the 3D uniform electron gas. 
Let us consider two adjacent electrons with the center-of-mass and relative coordinates, 
X = (x\ + X2}/2 and x = x\ — X2, respectively. Focusing on the x dependence the 
Schrodinger equation can be written as 



d 2 



dx' 



+ 



A(A-l) 



<p(x)$(X, x 3 , 



(E-H')<p(x)$(X,x 3 ,...),(B.l) 



where H' contains the remaining terms in the Hamiltonian. Because E — H' is non- 
singular as x approaches zero, it is unimportant for small x. To lowest order in x we 
therefore have <p(x) = x x + . . ., from which immediately follows n(x) ~ x 2X for the PD, 
see Equation (|4.1| ). This can be concluded for A ^ 1 also directly from the many-body 
wave function $ ~ U i<: jX^ [Q and for A = 1 from Equation (|4.11| ). 



A similar treatment of the asymptotic large k behavior of the momentum 
distribution n K seems to lead in Equation ( |4.3cj) to the conclusion 7 = 2A + 2. This 
corresponds to n K ^ c 



(7(0) /fir for the 3D uniform electron gas |40|, ff6| . [48 




Figure 1. Bulk energy e (solid), kinetic energy t (dashed), and potential energy v 
(dotted) plotted as functions of interaction strength parameter v. Thin lines denote 
the results of the Hartree-Fock approximation, thick lines are exact. The thin dashcd- 
dotted line indicates the "fall- into- the-origin" at v = 0. 
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Figure 2. Dimensionless PD g(y) = n{x\2)/n 2 as a function of the dimensionless 
inter-particle separation y — kpxi 2 for v = (dashed), 1/2 (solid), and 3/2 (dotted). 
The thin line is a guide to the eye only. 
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Figure 3. Fermionic momentum distributions n K vs. n = k/kp with v — (dashed), 
1/2 (solid), and 3/2 (dotted). Inset: n K for v = 3/2 and L — 401 computed with 
N = 21, 41, 81, 121, 161, 201, 241, 281, 321, 361, and 401. The data for N = 21(o), 
41(D), 81(o), and 121(x) do not show any density dependence whereas the larger 
density data (lines) do. 




Figure 4. Static structure factor S(q) = 1 — h(q)/n for v = (dashed), 1/2 (solid), 
and 3/2 (dotted). Inset: The three curves do not coincide at a single point close to 
q w 1.72. 




Figure 5. Particle-number fluctuation (ANx) 2 /Nx in domains X of the CS model 
after Equation for v = (dashed), 1/2 (solid), and 3/2 (dotted). The dashed- 
dotted line corresponds to (ATV^) 2 /Nx for strict correlation. || 
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Figure 6. Correlation measures based on (a) the bulk of the momentum distribution, 
(b) the tail of the momentum distribution, and (c) the PD shown as functions of v. 
The solid lines are guides to the eye only. The thin dashed-dotted line indicates the 
"fall-into-the-origin" at v — 0. 




Figure 7. Correlation entropy ( |5.6| ) for Fermions (solid line) as a function of density 
at v = 3/2. The dashed-dotted line corresponds to s obtained for the discrete CS 
model. 
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Figure 8. Kinetic energy t as computed from the Hellman-Feynman theorem (2.5) 
(dashed line), and t from Equation ( |2.6|) (solid lines) for Fermions (thick line) and 
Bosons (thin line) as a function of density at v = 3/2. 
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Figure 9. Bosonic momentum distributions n K vs. k — k/kp with v = (dashed), 
1/2 (solid), and 3/2 (dotted). Inset: n K for v = 3/2, L = 401 and particle numbers 
identical to the inset of Figure |. The data for N = 21 (o) , 41(D), 81 (o), 121(x), 
161 (+), and 201(*) do not show any density-dependence effects. 



